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ABSTRACT 



The mid-infrared properties of pre-planetary disks are sensitive to the temperature 
and flaring profiles of disks for the regions where planet formation is expected to occur. 
In order to constrain theories of planet formation, we have carried out a mid-infrared 
(A = 10.7/xm) size survey of young stellar objects using the segmented Keck telescope in 
a novel configuration. We introduced a customized pattern of tilts to individual mirror 
segments to allow efficient sparse-aperture interferometry, allowing full aperture syn- 
thesis imaging with higher calibration precision than traditional imaging. In contrast 
to previous surveys on smaller telescopes and with poorer calibration precision, we find 
most objects in our sample are partially resolved. Here we present the main observa- 
tional results of our survey of 5 embedded massive protostars, 25 Herbig Ae/Be stars, 
3 T Tauri stars, 1 FU Ori system, and 5 emission-line objects of uncertain classifica- 
tion. The observed mid-infrared sizes do not obey the size - luminosity relation found 
at near-infrared wavelengths and a companion paper will provide further modelling 
analysis of this sample. In addition, we report imaging results for a few of the most 
resolved objects, including complex emission around embedded massive protostars, the 
photoevaporating circumbinary disk around MWC 361A, and the subarcsecond binaries 
T Tau, FU Ori and MWC 1080. 

Subject headings: accretion disks — radiative transfer — instrumentation: interferom- 
eters — circumstellar matter — stars: pre-main sequence — stars: formation 



1. Introduction 

Theories of planet formation rely on observational estimates of preplanetary disk initial con- 
ditions in order to make predictions. By connecting detailed and accurate disk initial conditions 
with the observed diversity of exoplanetary systems, these theories promise to explain how, when 
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and where exoplanets form under a wide range of circumstances. Unfortunately, direct measure- 
ments of the initial conditions are very uncertain since these regions remain largely unresolved 
by conventional observing techniques. Indeed, disk theory is still undergoing active revision and 
new measurements more often reveal puzzling surprises rather than build confidence in pre-existing 
models. 

Pioneerin g near-infrared observations of y oung stellar ob.iect s (YSOs) using long-baseline in- 



terferometry flMillan-Gabet et all Il999l . 12001 



Tuthill et al. 



200 ll ) found disk emission sizes that 
were much larger than expected from the conventional disk theory. This has been successfully 
interpreted in terms of an optically-thin cavity surrounding the s tar, in contrast to previous as- 
sumptions of an optically-thick disk e xtending;; to stellar surface (IMonnier fc Millan-Gabetl 120021 : 
Natta et al.ll200ll : iMonnier et al.ll2005l ). While most work has been done for the intermediate mass 
Herbi g; Ae/Be stars, this g:eneral pattern is seen even for the young sun analogues, the T Tauri 
stars dAkeson et allbood : Isisner et al.|[2005l: iMillan-Gabet et allboOTl V 



The introduction of this large inner cavity to explain the near-infrared sizes radically changed 
the inner disk structure in models of these systems, requiring a "puffed-up inner rim" and even 
a shadowed region behind the rim in some models (j Pullemond et al. 2001 ). In order to probe 
the disk structure and surface temperatures at these larger radii, we need to observe at longer 
wavelengths. Mid-infrared size measurements using interferometric techniques promise to probe 
the terrestrial-planet forming zones around young stars, characterized by temperatures around 
300K. Until recently, most attempts to re solve YSO targets in the mid-IR have failed due to lack 
of angular resolution (e.g.. lLiu et al.lll996l ). 



Hinz and collaborators have re solved a few targets using a novel null ing interfometer BLINC 



Liu et al. 



20051 . 120071 ) . although most targets 



on a single 6.5m-aperture telescope rtHinz et alJ I I2OO1I : 
were unresolved. Clearly, longer baselines will be needed to increase the angular resolution; the ISI 
interferom eter was the first l ong-baseline interferometer to resolve thermal emission from a YSO 
(LkhalOl; iTuthih et al]l2002l l. 



A radical technological advance was recently made by the MIDI i nstrument on the VLTI 



interferometer, measuring seven mid-infrared sizes of Herbig Ae/Be stars (Leinert et al. 20041 ) and 



even taking a spatially-resolved spectra across the silicate feature (jvan Boekel et al.ll2004l ). With 
baselin es ranging from 1 0m to 200m, the VLTI-MIDI system (and also the Keck Interferometer 
Nuller; IStark et al.ll2009l ) is poised to revolutionize studies of YSO disks in the mid-infrared. 



In this article, we describe an experiment employing one of the world's largest single-aperture 
telescope, the Keck-1 telescope. We have optimized the calibration precision of our survey by pro- 
gramming the individual mirror segments of the Keck to create multiple non-redundant interfero- 
metric arrays. In this method, we allow precise calibration of changing atmospheric conditions while 
instantaneously sampling full 2-dimensional sky angles, allowing detection of disk elongations and 
asymmetries with unprecedented precision. Furthermore, we can also survey sources very quickly, 
allowing a large sample to be observed in only a few nights. Our study complements longer-baseline 
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observations since "short" baselines are needed to distinguish between la r ge-scale disk emission an d 
emission from the compact "puffed- up inner rim" (Ivan Boekel et al.ll2005l : iTannirkulam et al.ll2007l ) . 



Here we report the full results of our diffraction-limited size survey of bright young stellar 
objects, finding most of them to be resolved and a few to be elongated. In addition to our imaging 
results, we can investigate the mid-infrared size - luminosity relations for YSOs in detail for the 
first time. Unfortunately, the recent decommissioning of the only mid-infrared camera on the Keck 
telescope marks the end of this short-lived but exciting experiment. 



2. Observations 

2.1. Description of segment-tilting experiment 

While mid-infrared imaging with 10-m class telescopes is often considered "diffraction-limited," 
precise measurements of partially-resolved objects are difficult due to changes in seeing between 
observing a target object and a Point Spread Function (PSF) calibrator. In order to circumvent 
these difficulties, we designed a novel interferometer experiment whereby individual segments of 
the 36-segment Keck telescope were reconfigured to form four interferometric arrays. Since each 
segment (1.8m across) is much smaller than the typical coherence length of the atmosphere (rg =5- 
10m at A =10^m) the results are more robust to variations in seeing conditions. The general 
approach of apodizing a telescope for interferometry is known as "aperture masking" and the 
technique has well-established advantages over speckle interferom etry and adaptive o ptics for bright 

and the more 



targets. We refer the interested reader to the pioneeri ng work of 



recent implementation on Keck for further discussion (iTuthill et al 



Haniff et al. 

booni ). 



In this work, we formed five interferometric sub-arrays by tilting the Keck primary mirror 
segments in a customized pattern. Figure [1^ shows the Keck segment map color-coded by the 
applied tilt in a pattern known as "pattern 6" (while other tilt patterns were available, only pattern 
6 was used in this work). These tilts caused the stellar image to split into 5 separate images as 
illustrated in Figure [lb. Four of the images (each with light from 6 segments) were arranged in a 5" 
X 5" square which fit perfectly on the 10" x 10" field-of-view of the Long; Waveleng;th Spectrometer 



(LWS) camera on Keck-1, a 128x128 array with 0.83^^ platescale (ICampbell fc Jonesll2004l : IPerrin 



20061 ) . The four images were used for science and were measured during data acquisition. The 
segments contributing to each of the four "science" arrays were arranged in a non-redundant pattern 
so that the fringe visibilities of all baseline pairs could be extr acted using; Fourie r analysis (see 
discussion of redundancy in aperture masking interferometry in lTuthill et ahlBoOol ). Because not 
all segments could be incorporated into four non-redundant patterns, a fifth image was created 
using light from the 12 unused segments and this light was projected 15^^ away from the detector. 

This approach to aperture masking has a few significant advantages over those employing 
opaque masks placed in the telescope pupil plane (as done for all previous aperture masking work). 
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Firstly, since light from 24 of the 36 segments are used for science, the cohecting area of the Keck 
primary is efficiently used (compare to near-infrared Keck aperture masking which passes 1-10% 
of the incident flux). Secondly, each of the four sub-arrays generally has different (u,v) Fourier 
coverage, which is essential for high fidelity aperture syn thesis imafflng. Indeed, we have alread y 



presented hig!;h-re solution images of dus ty evolved stars (IWeiner et al.ll2006l: llr eland et alJ 120071 ) 



Wolf-Rayet stars flRajagopal et all 120071 ). and one YSO (IMonnier et al.ll2008l ) with unprecedented 
angular resolution. Thirdly, tilting segments is operationally straightforward and does not require 
any hardware modifications - we expect this technique to find application for the next generation 
of segmented telescopes, e.g.. Gran Telescopio Canarias, Cornell Caltech Atacama Telescope, and 
the Thirty Meter Telescope. 



2.2. Segment-tilting procedure 

The segments were tilted by sending offsets to the 108 actuators which support and control the 
Keck primary mirror. For all segments in a subarray, offsets were calculated so that incoming light 
was reflected to the desired off-axis location in the telescope image plane. Note that in addition 
to the tilts, a large piston offset was needed to be applied to the segments in order to preserve 
phase coherence at each of the four new pointing origins. In essence, this extra piston term was 
necessary so that each subset of segments conformed to a new parabolic surface. These piston 
offsets can amount to ^100/xm, critical to correct so that interferometric fringes fall within the 
coherence envelope of the light (typical narrow-band filters have ^ ^ 7). 

In normal Keck operation, offsets between mirror segments are continuously monitored and 
maintained by the Keck Active Control System (ACS), employing capacitive edge sensors. However, 
our experiment required the edge sensors to operate with a lower gain setting due to the much larger 
gaps between segment edges encountered here. Unfortunately, the sensor response was not well- 
calibrated under these conditions, and so the actual tilts and pistons realized in practice only 
approximate the desired configuration (typical errors were ^" ). Fortunately, the introduced errors 
were found to be mostly corrected by employing so-called "focus mode" to the primary. "Focus 
mode" approximates telescope focus by changing the overall curvature of the Keck primary mirror 
surface. It is a well-known characteristic of the Keck ACS system that changing sensor gain settings 
introduce focus-mode error (Richard Cohen, private communication). 

After correcting focus- mode errors, there were still residual tilt mismatches that affected some 
segments greatly. Clearly, if the light patterns from individual segment do not overlap well, the 
fringe power will be reduced and become sensitive to small changes of mirror figure (say, as a 
function of elevation). In order to optimize beam overlap, we applied a final correction using 
"optical feedback." That is, we took a set of data at the beginning of each night using our best 
focus-mode alignment and analyzed the resulting fringe patterns to estimate correction terms for 
each segment. The algorithm we used employed phase slopes in the Fourier transform of the speckle 
patterns and would not have worked in the case of a redundant array of segments. We then applied 
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these small corrections as a perturbation to the existing ACS offsets. This entire procedure took 
approximately 10 minutes and resulted in excellent alignment, with individual segments overlapping 
with j^^^ precision, approximately 0.1 -^-^ — -. See Figure El showing short-exposure fringe images 
before and after the final "optical tweaking." 

It was found that the segment alignment could change by up to ^0.5" as a function of telescope 
elevation and other operational variables (e.g., temperature), necessitating re-calibration. In gen- 
eral, we solved this by having separate actuator calibration "snapshots" as a function of elevation 
and using calibrator stars near to each target in time and elevation angle. For bright objects one 
can use post-processing to measure the tilts of each segment after the fact (akin to the "optical 
tweak" algorithm discussed above) and apply correction factors. We did not use this extra cali- 
bration step in this paper since the correction is signal-to-noise dependent and many of the YSOs 
were too faint for robust implementation. 



2.3. Data collection methodology 

We employed standard chop-nod observing methods to minimize noise from fluctuating back- 
ground. The secondary mirror was chopped at 5 Hz with a 10^^ throw along the telescope azimuth 
direction. In general, each saved data frame had a total integration time of 90 ms and we col- 
lected approximately 350 individual chop-nod sets. We sometimes employed longer integrations 
and slower (2.5 Hz) chop cycles for faint targets, and also collected more chop-nod sets to im- 
prove signal-to-noise ratio. The 90 ms integration times were sufficient to essentially "freeze" any 
fringe blurring due to atmospheric turbulence at these wavelengths, although telescope wind-shake 
can still cause calibration problems; fortunately, wind-shake was minimal during the observations 
presented here. Each target was observed two independent times, bracketed by observations of 
point-source calibrators located nearby in the sky. Table 1 contains the complete observing log, 
including the target names, dates, times, and calibrators used. We observed all of our targets in 
only one filter, the 10.7/xm filter (passband 9.92-11.47^m). 



2.4. Data reduction 



The bulk of the data analysis processing was carried out by the same aperture masking code 
developed for the near -infr ared version o f the K eck masking experiment, which is best described in 
Monnier et al.l (Il999al ) and lTuthill et al.l (120001 ). The only major modification to the data pipeline 
was a change meant to address significant camera noise in the power spectra of our short exposure 
images. Here we used the nod sequences (which did not have any star light) to estimate a con- 
temporaneous bias power spectrum. This bias subtraction proved necessary and key to extracting 
reliable results for the faintest targets in our sample. 



The overall observing procedure was introduced in the first segment-tilting publications (IWeiner et al 
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20061 : iRajagopal et alJ 120071 : iMonnier et al.ll2008l ). Essentially, each subarray pattern (field-of-view 
3^' ) was analyzed using Fourier techniques, resulting in fringe visibilities for each baseline and clo- 
sure phases for all triangles. These values were calibrated using point-source reference stars. Some 
objects are very resolved and the shortest inter-segment spacing is approximately 2 meters long, 
meaning we are missing a lot of short-baselines. In order to make up for this, we also analyzed the 
same data frames using baselines shorter than 1.8 meters, which correspond to sampling the many 
redundant baselines within a segment. These short baselines are useful when reconstructing a wide 
binary companion during apertures synthesis imaging and helps to constrain information on the 
percentage of the light coming from a large-scale "halo" for some targets. 

The measurements from the four simultaneous (inter-segment) patterns and the results of 
the "short-baseline" (intra-segment) measurements are merged together and used in subsequent 
analysis. We used observations of the known binaries MWC 1080 and SVS 20 to validate our 
position angles and pla te scale. The ca librated data are saved in the OI-FITS data format for 
optical interferometry (iPauls et al.ll2005l ) and are available upon request. 



2.5. Analysis of Systematic Errors 

Since we expected some of our target sample to be "unresolved" it was critical to establish 
reliable upper limits in those cases. This required an ex tensive analysis of sy stematic errors, as was 



recently done for near-infrared aperture masking data (IMonnier et al.l 120071 ). 



The best method for reliable error estimation is to obtain multiple independent datasets for each 
object. This requirement was a driving consideration in our observing strategy and we obtained 
multiple observations for all our targets as can be seen in Table 1. By analyzing the consistency 
between the multiple independent measurements we can provide robust error analysis. 

Figure [3] Shows a consistency plot for all the data of the Keck segment tilting experiment using 
Pattern 6 - this includes many objects not part of this YSO survey but still useful for the analysis 
of systematic errors. Here we graph the size measured for an object at one time against the size 
measured for the same target at another time (and nearly always using an independent calibrator 
sequence). To create this graph we removed known binaries (MWC 1080, T Tau, FU Ori) and 
the error bars were estimated by analyzing the 4-subarrays of pattern 6 using bootstrap sampling 
techniques. Based on the variation in measured diameter as function of mean diameter, we can 
quantitatively estimate our 2-sigma detection limit as FWHM 35 milliarcseconds - our confidence 
limits as a function of mean diameter are also included in Figure [3l 



As was found in IMonnier et al.l (120071 ). the calibration precision is non-linear in the sense that 
the fractional error is very large for small sizes while being greatly reduced for larger sizes. Our 
final fitting results (both 1-D and 2-D) used bootstrap sampling of the different epochs and sub- 
arrays to capture calibration and systematic errors. In the case that two independent measurements 
disagreed dramatically, we removed this target from our sample. As can be seen in Figure [HI this 
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rarely occurred and did not significantly impact our sample size. 



2.6. Aperture synthesis imaging 



Aperture synthesis imaging can be carried out for the most resolved targets, useful for de- 
tecting faint diffuse emission and for det ecting; binary companions. We used the publicly available 



BSMEM image reconstruction software (lBuscherlll994l : iLawson et al 



200' 



-zi. .20061 ) for aperture syn- 
thesis imaging. This program was discussed recently i n iMonnier et al. | (20081) and wa s validated 
in more detail against the MACIM fllreland et Mlbflod) algorithm in|zhao et all fl2008l ). BSMEM 
uses the maximum entropy met hod (Gull fc SkillingI 1983 : Narayan fc Nityananda 1986 ) and is sim- 
ilar t o the VLBMEM prog r am (ISivial 119871) extensiv ely used in the near-IR Keck masking project 
(e.g., iMonnier et al.lll999bl : iTuthill et al.lll999l . 12000). While most young stellar objects were only 
partially resolved here, a few sources were large enough (on the sky) to be imaged and these are 
discussed in 



3. Results 



Table 2 contains the complete source list for this work, including 5 embedded massive pro- 
tostars, 25 Herbig Ae/Be stars, 3 T Tauri stars, 1 FU Ori systems, and 5 emission-line objects 
of uncertain classification. We have also included other important characteristics of our sample, 
including V, J, H, and K magnitudes, our measured lO.T^m flux within aperture of 3^^ diameter, the 
IRAS 12 and 25/xm flux, and notation regarding the presence of a close binary companion (based 
on the literature). Subsequent tables and figures maintain the groupings and target orderings as 
presented in Table 2. 



3.1. Visibilities 



Figure 4 shows the squared- visibility data for all our targets as a function of baseline in meters. 
In this figure, we have performed (u,v) averaging on scales of 0.45 meters. This averaging is done 
in the (u,v) plane so baselines are only merged when sharing similar lengths and position angles. 
Here we did not average together results from the different independent observations, but rather 
overplotted each separate epoch detailed in Table 1. Thus, the observed scatter represents the 
combination of observing errors plus any position-angle dependency of the source structure. For 
instance consider the visibilities of v892 Tau and MWC 349A which vary for the same baseline 
length outside the error bars - indeed, these two objects are extremely asymmetric. Also consider 
that MWC 1080 and T Tau are both well-known sub-arcsecond binary stars - this is obvious in 
the scattered visibility data. 
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Figure 5 shows the squared-visibihty as a 2-dimensional function of (u,v) coordinates, a useful 
presentation format for detecting elongations. For instance, the extreme asymmetry for MWC 349A 
and v892 Tau, inferred from the 1-D visibility curves, is quite evident here. Also, the large "halo" 
structure around MWC 361A, indicated by the sharp drop of visibilities as short baselines and then 
a plateau, can be easily spotted. The binary signatures in T Tau, FU Ori, and MWC 1080 are 
evident too, although it is difficult to interpret since the sinusoidal patterns are under-sampled, 
requiring analysis through imaging or model-fitting to extract the binary component separation. 



3.2. Fitting Results 

We have characterized the emission sizes of all targets by fitting the visibility data with both 
a circularly-symmetric (1-D) Gaussian model and an elliptical (2-D) Gaussian model. In both 
cases, we also allowed for a large-sc ale "halo" componen t, similar to the model used for near- 



infra red fitting of YSO visibilities in iMonnier et al.l (|2006l ) which was motivated by earlier work 



(e.g. iLeinert et al.l 12001 



Table 3 contains the complete fitting results from our work. Error bars were estimated through 
bootstrap sampling of the multi-epoch data as well as from the different sub-arrays of "Pattern 6." 
For each fit (1-D and 2-D) we also report the best fit xi normalized by degrees of freedom. Thus, 
if the Xu value is greater than ~1 then we expect this model to be an incomplete description of the 
data. For instance for MWC 349A, the xt is 4.0 for the circular-symmetric fit but decreases to 1.1 
for a 2-dimensional Gaussian - we already discussed above that this object was clearly asymmetric. 
In most cases when the Xu is ^1 for the 1-D fit, then the 2-D fit is only marginally better and 
that there is no statistically-significant ellipticity detected. Some targets have high xl even for 
the 2-dimensional Gaussian model, indicating a more complicated object. This is true for all the 
embedded YSOs which are generally the most resolved - these objects are ideal targets for aperture 
synthesis imaging in the next section §3.31 



3.3. Imaging 

While most of the analysis in this work and the companion paper relies on the size-fitting 
results, we also report imaging results for the most resolved targets. We used the BSMEM image 
reconstruction program (see §2.6^ on all targets and inspected the results. In Figure El we collected 
the images of all the objects that showed structures beyond a simple, partially-resolved Gaussian 
profiles and/or were larger than 70 mas in size. Note there are residual imaging artifacts at the 
level of a few percent of the peak, especially evident for the binary sources MWC 1080, T Tau and 
FU Ori. Notes on individual objects follow. 

AFGL 490 has extension to south. 
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Mon R2 IRS 3 shows nebulosity along a N-NE to S-SW axis a nd a companion a t p — 836 mas, 
9 — 16.7°, similar position to the companion first reported by iMcCarthyl (Il982l ) of separation 
870 mas at PA 13.5°. 

AFGL 2136 has a slight extension to west-southwest that appears in two independent sets. 
AFGL 2591 emission has slight extension to west. 
S140 IRSl shows significant nebulosity to the South. 

AB Aur is a prot otypical Herbig Ae star which has been claimed to be elongated nearly 2- 
to-1 along PA 30° bv iLiu et all (|2005l ). We rule out this large level of elongation (supported by 
4 separate observations over 3 different observing nights), although we do find a slight extension 
along this position angle in our image. Overall the emission is circula rly symmetric, consis tent with 
a nearly face-on geometry and in agreement with the conclusions of iMarinas et al.l (j2006l ). 



v892 Tau was resolved into a circumbinary disk and was discussed in iMonnier et all (|2008l ) . 
Here, we reproduce the image. 

LkHa 101 shows a s light elongation alo ng North-South direction, similar to reported asymme- 
try in the near-infrared (jTuthill et al.ll2002l ). Note we do not have fine enough angular resolution 
to resolve the dust-free cavity in the center reported by these authors. 

R Mon is resolved approximately as much as AB Aur, showing mostly symmetrical emission. 



Z CMa is a well-known close binary (^100 mas) recently imaged by iMillan-Gabet fc Monnier 

(120021 ) . Here, the mid-infrared emission appears extended towards position angle ^128°, which 



is the direction of the FU-Ori type companion. This extension could represent emission from a 
circumbinary disk or direct IR emission for the FU Orionis component to this system. Most likely, 
the bulk of the mid-IR emission is from the Herbig Be component of the system and not the FU 
Ori object. 

MWC 349A is very elongated and the imaging sh ows a symmetric str ucture oriented along PA 



95°, very similar to the PA 100° reported in the NIR flDanchi et all 12001 



MWC 361A was found to have a large "halo" containing 45% of the flux at 10.7^m. This 
halo was over-resolved on our shortest baselines. Our shortest baseline data suggested the halo 
was elongated north-south on arcsecond scales and this is confirmed by independent (standard 
full-aperture) LWS images taken at 11.6/xm (SiC filter) and 17.65/xm. These image s are shown i n 
Figure El a full description of these observations and the data reduction appears in IPerrinI (2006). 
Note the N-S elongation of the extended emission matches the orientation of the binary orbit of 
MWC 361 A measured by iMonnier et al. I (l2ond ). strongly suggesting this halo is the remnant of a 
circumbinary disk. It is remarkable that the mid-IR disk emission is >20x bigger than the semi- 
major axis of the binary, indicating a huge disk gap in this system much larger than can be produced 
by disk clearing through dynamical interactions with the inner binary. Given the early spectral 
type of the primary star in this system, we suggest that phot oevap oration of the circumbinary disk 
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is well underway. 

MWC 1080 is binary with p = 764 mas, 9 

fll997l ) of p = 760 mas, 9 = -93°. 



-91.1°, comparable to the result of iLeinert et al 



The T Tau companion to the south (which is itself a close binary) is detected at p = 639 mas, 
9 — —170.9° We attempted a more sophisticated analysis to look for evidence of the tertiary to 
this system. Indeed , we find that that the visibility and closure phase data can not be fit with a 



simple binary star (ISkemer et al.l 12008 ^). However, there are calibration issues regarding field-of- 



view limitations caused by bandwidth-smearing that we have not corrected here. This issue will be 
tackled in a future paper. 

We detec t a companion to F U Ori at separation 488 mas and PA of 163.3° E of N, confirming 
the report by I Wang et al.l (120041 ) of a potential companion at separation 500 mas at PA 161°. As 
for T Tau, we can not report a reliable fiux ratio due to unresolved issues with bandwidth smearing 
effects that have not yet been corrected. 



3.4. Comparison to Literature 



Our resul ts are g!;enerally i n agreement w i th res ults from the BLINC experiment by Hinz and 
collaborators flHinz et al.ll200ll : lLiu et al.ll2005l . 120071 . excepting the elongation of AB Aur), although 
we note most of these workers' measurements were upper-limits. We also comp ared our results t o 
the long-baseline interferometry VLTI data of iLeinert et al.l (120041 ) and others ( Quanz et al.ll2006l ). 



finding that the Gaussian fits do not always agree between the two methods (our sizes typically 
larger). The size discrepancies can be understood if the long-baseline VLTI measurements are more 
sensitive to the properties of the hot inner wall while th e short-baseline (single) Keck data can only 



prob e the large-scale, outer fiare d disk (see discussion in Ivan Boekel et al.ll2005l : iTannirkulam et al 



20071 : iMillan-Gabet et all 120071 ). These observed differences appear to confirm the expectations 
from models that mid-infrared emission comes from (at least) two distinct spatial scales and we will 
need measurement from a wide range of scales in order to reconstruct the actual disk temperature 
profile. 



4. Analysis 

4.1. Size — Luminosity Diagram 

A very tight correlation has been found between the near-infrared siz e of a young stellar object 
disk a nd the luminosity of the host star for Herbig A e and late Be stars (IMonnier fc Millan-Gabet 



2002 



Monnier et al. 



20051 : iMillan-Gabet et al.l 120071 ). This simple correlation arises because the 



NIR emission only comes from the hottest dust near the evaporation front and the location of the 
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evaporation front is primarily sensitive to a single observable, the central luminosity. 



The mid-infrared emission from young stars is expected to come from up to three regions: 
the hot inner wall loc a ted at the dust evaporat ion radius, a thin surface layer on the disk (e.g. 
van Boekel et al.ll2005l : iTannirkulam et al.ll2008l ). and perhaps a circumstellar envelope ("halo"). 
Thus, the mid-infrared size and flux density will depend on multiple factors in addition to the 
central luminosity, depending most sensitively on the dust size distributions both vertically (setting 
temperature) and radially (controlling disk flaring). While it is beyond the scope of this paper to 
investigate these effects, we will discuss them generally in the context of the mid-IR size - luminosity 
diagram. 

Figure [8] shows the mid-infrared siz e - luminosity diag:ram fo r the stars in our sample, con- 
structed using the same procedure as in iMonnier fc Millan-Gabetl (120021 ) . Our sample represents 
a factor of 3 times increase in sample size over prev ious work and spans a l a rge r ange of YSO 
luminosity. We used distances and luminosities from lAcke fc van den Anckerl (120041 ) when avail- 
able. In this diagram, we flnd a crude correlation of size with luminosity showing disk emission 
with characteristic temperatures of 250-900 K spanning about 5 orders of magnitude in luminosity. 
While the previously-published near-infrared size - luminosity diagram shows only about a factor 
of 2 scatter at a given luminosity, this mid-infrared relation shows much larger scatter, a factor of 
5 or more. The origin (s) of the increased scatter and the breakdown in the size-luminosity relation 
will be the subject of a companion paper (Tannirkulam et al., in preparation). Here, we will only 
make general comments. 



In ITannirkulam et al.l (120081 ). the large mid-infrared size difference between two "twin" Herbig 
Ae stars (AB Aur and MWC 275) was explained through the influence of small grains in the AB Aur 
disk at radii of ^7 AU. These grains absorb stellar radiation effectively, heating up the disk and 
causing flaring. The origin of these small grains could be due to collisions between planetesimals or 
the release of small grains as ice-cemented dust agglomerates break apart at about ~10 AU. This 
example illustrates how variations between targets in the radial dust size distributions can explain 
some of the scatter in this diagram. Other modiflcations to the dust dis tribution, such as dust 
grow t h and settling, can als o cause signiflcant changes in mid-infrared SED (iDullemond fc Dominik 
2004l : lD'Alessio et al.ll2006l ). 



Another source of scatter is the presence of binaries. We have indicated close binaries in 
Tableland in Figur e [8l For instance, th e emission from v892 Tau was recently imaged to be from 
a circumbinary disk (IMonnier et al.ll2008l ) and we see in Figure [8] that v892 Tau is indeed oversized 
for its luminosity. 

We also note that the embedded YSOs all cluster near the top-right of the diagram at high 
luminosities and low temperatures. This is expected since the high optical depths toward these 
Class I objects mean their infalling envelopes are still optically-thick and we do not expect to be 
able to see into the warmer inner regions surrounding these protostars. Note that these regions all 
show complicated, often bipolar, structures (see Figure [6j) and that the sizes in this diagram refer 
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only to the central core of emission. 



4.2. Size vs IRAS color 



Leinert et alJ (120041 ) found some evidence that the physical emission scale for the mid-infrared 
emission (in Herbig Ae stars) was correlated with the IRAS 12-25 micron color, defined as - 
2.5 log Fiy(12/xm) / Fiy(25/xm). This is sensible since redder colors indicate cooler dust temperatures 
which naturally arise farther from the central star. Unfortunately, this diagram is difficult to use 
when your target sample spans a large luminosity range. This is because the radius for dust at 
given temperature depends on the root of central luminosity (see size-luminosity diagram above). 
Neither size-l uminosity nor t he size-color diagram can simultaneously take into account both of 



these effects. iLiu et al.l (120071 ) also explored mid-IR size differences as a function of the some SED 



parameters but his work also suffers from these same limitations. 

Despite these drawbacks, we have produced a similar diagram as a way to present our basic 
results. Figure [Qfei shows our measured Gaussian FWHM size (in milharcseconds) vs. the IRAS 
12-25 color. We chose to use the distance-independent angular size instead of the physical size 
here, just so that our quantities did not depend on the uncertain distance estimates. Within each 
grouping, we do see some correlations. The Herbig Ae stars do show a correlation, with the striking 
result that nearly all the objects with colors bluer than 0.25 were unresolve d by our survey. Th e 



embedded objects also might show a weak correlation. Just as the diagram in lLeinert et al.l (120041 ) . 
this size-color diagram does not conserve a target's position if you take the same disk temperature 
profile and change the luminosity. We only present this data to illustrate some trends within the 
sample, but caution against further use of this sort of diagram in subsequent work. 

In order to overcome the luminosity and distance-dependencies for this sort of diagram, we 
introduce a new diagram in Figure [9)3. Here plot the luminosity-normalized ring radius versus 
iras 12-25 color. This is defined as the physical ring radius (AU) divided by the square-root of 
luminosity (in solar luminosities). This quantity is both independent of distance and also naturally 
accounts for the expected scaling of emission scale with luminosity (for similar disk temperature 
profiles). Again, we see a significant correlation for the Herbig Ae stars, although other objects 
appear to scatter in this diagram. The location of each YSO system in this diagram is diagnostic 
of the disk flaring and temperature profiles and this diagram will prove useful for future work. 

In the next paper in this series, we will be carrying out a radiative transfer study of YSO 
disks in order to understand the scatter in the size-luminosity diagram and how the deviations 
correlate with observables such as stellar luminosity and the [12]- [25] micron color. This new work 
will motivate improved diagrams to use for plotting basic observational data and will explore the 
potential for using luminosity-normalized sizes (radius AU / sqrt Luminosity) and surface-brightness 
relations. 
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5. Summary 

We have presented a large survey of mid-infrared sizes of young stellar objects using the Keck 
Telescope. We can reliably resolve disks with size scales down to 35 mas (5 AU at 140pc), roughly 
8x smaller than the formal diffraction limit. We have overcome the traditional calibration problems 
that have plagued previous surveys by reconfiguring the Keck segments into multiple non-redundant 
interferometric arrays. 

We have presented a detailed description of the experimental methodology and validated our 
calibration precision. We have measured the characteristic sizes of all our targets using a Gaussian + 
Halo model and have presented aperture synthesis imaging when possible. Our atlas of diffraction- 
limited imaging has revealed the complicated emissions around a sample of Class-I embedded YSOs 
and also discovered the remarkable photoevaporating circumbinary disk around MWC 361A. 

Lastly we have presented the most complete size - luminosity diagram for YSOs in the mid- 
infrared. Notably we found the mid-IR size - luminosity relation shows a factor of 5-10 scatter 
for a given luminosity much larger than the tight correlation seen in the near-infrared. This large 
scatter can partially be understood in terms of the influence of small dust grains that can be create d 



in collisions between planetesimals in these young planet-forming disks (iTannirkulam et al.ll2008l ). 
Disk clearing by binary companions and emission from remnant dust envelopes may also contribute 
to this scatter. With the help of a large grid of radiative transfer models, we will be fully exploring 
how the mid-infrared sizes correlate with disk properties in the next paper in the series. 

While this particular experiment can no longer be carried out due to the absence of mid-IR 
instrumentation on the single Keck telescopes, we hope our method can be applied on current and 
future proposed segmented telescopes, such as the Gran Telescopio Canarias, the Cornell Caltech 
Atacama Telecope, and the Thirty-Meter Telescope. 



We have appreciated the involvement of C.H. Townes in this project, especially for his contri- 
butions to the initial work. We give special thanks to Mark Kassis and Randy Campbell for su- 
perb observing assistance and helpful discussions concerning experimental details. We also thank 
Jayadev Rajagopal for his help at the telescope. We acknowledge partial support from NASA 
KPDA grant 1267021, NASA Origins grant NNG05GI80G, nsf-ast 0352728, and the NASA Michel- 
son Fellowship program (MJI). This research has made use of the SIMBAD database, operated at 
CDS, Strasbourg, France. The data presented herein were obtained at the W.M. Keck Observatory, 
which is operated by Caltech, University of California and NASA. WMKO was made possible by 
the financial support of the W.M. Keck Foundation. The authors wish to recognize and acknowl- 
edge the very significant cultural role and reverence that the summit of Mauna Kea has always 
had within the indigenous Hawaiian community. We are most fortunate to have the opportunity 
to conduct observations from this mountain. 
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Table 1. Observing Log 



Target 


U.T. 

Date 


Time 


Galibrator Names^ 


51 Oph 


2005 


May 26 


10:05 


(7 Lib 




2005 


May 26 


10:20 


cr Lib, v3879 Sgr 


AB Aur 


2004 


Aug 30 


15:12 


i Aur 




2004 


Aug 30 


15:34 


i Aur 




2004 


Aug 31 


15:29 


i Aur 




2004 


Sep 01 


15:37 


i Aur 


AFGL 2136 


2005 


May 26 


12:47 


v3879 Sgr 




2005 


May 26 


13:06 


v3879 Sgr 


AFGL 2591 


2004 


Sep 01 


09:01 


X Oph 




2004 


Sep 01 


09:28 


X Oph 


AFGL 490 


2004 


Sep 01 


12:23 


a Tau 




2004 


Sep 01 


12:40 


a Gas 


FU Ori 


2005 


Feb 20 


07:28 


i Aur 




2005 


Feb 20 


08:33 


i Aur 


GW Ori 


2005 


Feb 20 


07:14 


i Aur 




2005 


Feb 20 


08:22 


a Tau, i Aur 


HD 142527 


2005 


May 26 


08:32 


GM Lup 




2005 


May 26 


09:12 


GM Lup 


HD 142666 


2005 


May 26 


09:49 


a Lib 




2005 


May 26 


10:11 


cr Lib 


HD 144432 


2005 


May 26 


08:56 


GM Lup 




2005 


May 26 


09:42 


a Lib 


HD 41511 


2005 


Feb 20 


05:24 


a GMa 




2005 


Feb 20 


09:21 


a GMa 


HD 45677 


2005 


Feb 19 


07:32 


a GMa, a Tau 




2005 


Feb 19 


07:32 


a GMa, a Tau 




2005 


Feb 19 


09:08 


a GMa 


HD 50138 


2005 


Feb 20 


05:50 


a GMa 




2005 


Feb 20 


09:16 


a GMa 


LkHa 101 


2004 


Aug 30 


14:45 


a Tau, i Aur 




2004 


Aug 31 


14:56 


a Tau, i Aur 




2004 


Sep 01 


15:33 


i Aur 




2005 


Feb 19 


07:00 


a Tau 


LkHa 234 


2005 


May 26 


14:54 


MO Gep 




2005 


May 26 


15:11 


MO Gep 


MWG 1080 


2004 


Sep 01 


13:25 


a Gas, MO Gep 




2004 


Sep 01 


13:28 


a Gas, MO Gep 


MWG 147 


2005 


Feb 20 


05:59 


a GMa, a Tau 




2005 


Feb 20 


08:54 


a GMa 


MWG 275 


2004 


Sep 01 


06:33 


v3879 Sgr 




2004 


Sep 01 


06:46 


v3879 Sgr 


MWG 297 


2004 


Aug 30 


09:24 


v3879 Sgr 




2004 


Aug 31 


07:39 


v3879 Sgr 


MWG 300 


2004 


Aug 31 


07:44 


v3879 Sgr 




2004 


Aug 31 


07:55 


v3879 Sgr 
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Table 1 — Continued 



Target 


U.T. 

Date 


Time 


Calibrator Names^ 


A/TTXT'/^ O /I O 


2UU4 Aug 61 


10:15 


XI Cyg 




zUU4 Aug dl 


10:43 


"VT r^^T^ -,^1 oon z^-,,^ 

XI Cyg, vl339 Cyg 




zUU4 Aug 61 


10:10 


XI Cyg 




zUL)4 Aug Jl 


10:53 


XI Cyg, vl339 Cyg 


MWU obl 


or>r>c i\/r, o/^ 

z(J(J5 May zb 


14:48 


MU Cep 




zUUo May zb 


15:03 


MU Uep 


A/TTXT'/^ /ion 


2UU4 Aug 61 


15:06 


i Aur 




ZL)LI4 Aug oi 


15:10 


i Aur 




ZLIU4 Aug Jl 


15:54 


a Tau, i Aur 




zL)U4 Aug 61 


15:58 


a Tau, i Aur 




zUU4 bep Ul 


15:06 


i Aur 


MWU Di4 


zUU4 Aug oi 


nn.no 

uy:Uz 


7 Aql 




zUU4 Aug oi 


nn. OQ 
Uy:zo 


7 Aql, R Lyr 




zUU4 bep Ui 


05:47 


7] bgr, vooTy bgr 




zL)L)4 bep Ui 


06:08 


Tj bgr, vo57y bgr 


Mon riz irib 6 


zUUo i^eb iy 


07:48 


OL CMa 




zUUo i^eb iy 


09:28 


OL CMa 


rli iau 


zUU4 Aug oi 


io:ib 


ex. Tau 




Onn/l A lire QI 

ZUU4 Aug oi 


io.io 


OL Tau 




zUU4 Aug oi 


15:49 


OL Tau 




zUU4 bep Ui 


15:11 


i Aur 


ID r^v A 
rt l^rA 


zUU4 Aug oi 


Ub: ib 


?7 Sgr 




zUU4 Aug oi 


06:34 


?7 Sgr 


R Mon 


zUUo i^eb zU 


05:30 


a. CMa 




zUUo reb zU 


05:33 


^ /^A /r^ 




onnc; t7^k on 
zUUO reb zU 


nn.nQ 
uy:Uo 


OL uivia 


Q1 /in TTD Q1 


onn/i Q^T^ ni 
zUU4 bep Ui 


1 1 . Qn 
ii:oU 


NiKJ Uep 




onn/i c^-.^ ni 
zUU4 bep Ui 


11:49 


MU Uep 


T Tau 


onnc T7^'U on 

zUUo i^eb zU 


06:39 


OL Tau, i Aur 




2005 Feb 20 


07:50 


a Tau, i Aur 


Z CMa 


2005 Feb 19 


07:56 


q; CMa 




2005 Feb 19 


09:33 


a CMa 


vl295 Aql 


2005 May 26 


14:04 


7 Aql 




2005 May 26 


14:19 


7 Aql 


vl685 Cyg 


2004 Aug 31 


09:32 


7 Aql 




2004 Aug 31 


10:58 


XI Cyg, vi339 Cyg 




2004 Aug 31 


11:02 


XI Cyg, V1339 Cyg 




2004 Sep 01 


09:52 


vl339 Cyg 


v376 Cas 


2004 Sep 01 


11:25 


a Cas 




2004 Sep 01 


11:44 


a Cas 


v645 Cyg 


2004 Aug 30 


11:28 


(3 And, 7 Aql, R Lyr 




2004 Aug 31 


10:38 


XI Cyg, V1339 Cyg 




2004 Aug 31 


11:19 


XI Cyg, V1339 Cyg 


v883 Ori 


2005 Feb 19 


07:38 


a Tau 




2005 Feb 19 


09:17 


a CMa 
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Table 1 — Continued 



Target U.T. Calibrator Names'^ 



Date 


Time 








v892 Tau 2004 Aug 31 


14:46 


OL 


Tau, 


i Aur 


2004 Aug 31 


14:51 


a 


Tau, 


i Aur 


2004 Aug 31 


15:34 


a 


Tau 




2004 Aug 31 


15:39 


a 


Tau 




2004 Sep 01 


14:58 


a 


Get 




2005 Feb 19 


07:07 


a 


Tau 





^All data taken using Pattern 6 (s ee JJT}. All calibrators were assumed to be as unresolved except for a Tau with UD 
diameter 19.8 mas (|Perrin et al.lll998h . 
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Table 2. Source List 



Target 


Alternate 


J2000 Coordinates^ 




J 


H 


K 


10.7/im 


IRAS12 


IRAS25 


Binary 




Name 


a 


5 










Jy^ 


Jyd 


Jy 


(Ref^ 


Embedded Young Stars (Class 


I Sources) 




















AFHT, 4Qn 




03 27 38.77 


+58 47 00.1 


- 


10.9 


8.1 


5.7 


58.3+ 5.8 


82.4 


278.0 


N 


IVKJIl JTLZ iJTLO O 




06 07 47.86 


-06 22 56.0 


- 


13.2 


9.8 


6.6 


109.9+11.0 


470.0 


4100.0 


Y(l) 


AFPT 91 "^fi 




18 22 26.38 


— 13 30 12.0 


- 


15.0 


12.7 


7.3 


37.1+ 3.7 


155.0 


574.0 


N 






zU zy 24.87 


+4U 11 19.4 


- 


14.3 


10.8 


6.6 


127.7+12.8 


439.0 


1110.0 


N 






22 19 18.28 


+63 18 45.8 


- 


12.3 


9.3 


6.1 


125.1+12.5 


308.0 


1540.0 


N 


Herbig Ae Stars 
























AJd Aur 


HD 31293 


04 55 45.83 


+30 33 04.4 


7.1 


5.9 


5.1 


4.2 


28.3+ 2.8 


27.2 


48.1 


N 


i\/r^Aro /ion 


HD 31648 


04 58 46.26 


+29 50 37.1 


7.7 


6.9 


6.3 


5.5 


15.4+ 1.5 


10.2 


10.3 


N 


riU i4zDOD 


vl026 Sco 


15 56 40.02 


-22 01 40.0 


8.8 


7.4 


6.7 


6.1 


5.9+ 0.6 


8.6 


11.2 


N 


rlJJ i4zoz / 




15 56 41.89 


-42 19 23.3 


8.3 


6.5 


5.7 


5.0 


11.1+ 1.1 


10.4 


21.2 


N 


TTTA -| A A AOO 

rlU i444oZ 




16 06 57.95 


-27 43 09.4 


8.2 


7.1 


6.5 


5.9 


9.8+ 1.0 


7.5 


9.4 


Y(2) 




HD 150193 


16 40 17.92 


-23 53 45.2 


8.9 


6.9 


6.2 


5.5 


18.9+ 3.8 


17.6 


18.1 


Y(3) 


51 Oph 


HD 158632 


17 31 24.97 


-23 57 45.3 


4.8 


4.9 


4.7 


4.3 


16.9+ 1.7 


15.7 


10.2 


N 


IVl W Aid 


HD 163296 


17 56 21.29 


-21 57 21.8 


6.9 


6.2 


5.5 


4.8 


15.3+ 1.5 


18.2 


21.0 


N 






19 01 53.68 


-36 57 08.2 


11.5 


6.9 


5.0 


2.9 


110.8+11.1 


111.0 


222.0 


N 


MWC 614 


HD 179218 


19 11 11.24 


+15 47 15.6 


7.2 


7.0 


6.6 


6.0 


20.6+ 4.0 


23.4 


43.6 


N 


vl295 Aql 


HD 190973 


20 03 02.51 


+05 44 16.7 


7.8 


7.2 


6.6 


5.9 


6.8+ 0.7 


7.2 


5.5 


N 


Herbig Be Stars 
























v376 Cas 




00 11 26.52 


+58 50 03.7 


15.6 


11.0 


8.4 


6.3 


21.1+ 2.1 


33.0 


88.1 


N 


v892 Tau 


Elias 1 


04 18 40.62 


+28 19 15.5 


15.3 


8.7 


7.0 


5.8 


34.1+ 3.9 


33.2 


99.9 


Y(4) 


LkHa 101 




04 30 14.44 


+35 16 24.0 


15.7 


8.3 


5.6 


3.0 


289.7+29.0 


362.0 


340.0 


Y(5) 


v883 Ori 


IRAS 05358-0704 05 38 18.10 


-07 02 25.9 


14.4 


9.3 


6.8 


5.2 


39.8+ 6.0 


52.5 


127.0 


N 


MWC 147 


HD 259431 


06 33 05.20 


+10 19 19.9 


8.8 


7.5 


6.7 


5.7 


10.0+ 1.0 


12.5 


20.2 


Y(6) 


R Mon 




06 39 09.95 


+08 44 09.7 


10.4 


9.7 


8.0 


6.4 


42.8+ 4.3 


54.7 


132.0 


Y(7) 


Z CMa 




07 03 43.17 


-11 33 06.3 


9.9 


6.5 


5.2 


3.8 


144.7+14.5 


127.0 


221.0 


Y(8) 


MWC 297 


NZ Ser 


18 27 39.53 


-03 49 52.1 


12.3 


6.1 


4.4 


3.0 


84.7+ 8.5 


159.0 


224.0 


N 


vl685 Cyg 


MWC 340 


20 20 28.25 


+41 21 51.5 


10.7 


7.9 


6.8 


5.8 


4.6+ 0.8 


- 


- 


Y(9) 


MWC 349 


vl478 Cyg 


20 32 45.53 


+40 39 36.6 


13.2 


6.2 


4.8 


3.2 


153.9+37.9 


179.0 


112.0 


Y(10) 


MWC 361 


HD 200775 


21 01 36.91 


+68 09 47.7 


7.4 


6.1 


5.5 


4.7 


8.2+ 0.8 


26.7 


76.8 


Y(ll) 


v645 Cyg 


AFGL 2789 


21 39 58.25 


+50 14 20.9 


13.6 


10.9 


9.2 


6.8 


57.4+12.0 


114.0 


219.0 


N 


LkHa 234 




21 43 06.82 


+66 06 54.2 


11.9 


9.5 


8.2 


7.1 


3.9+ 0.4 


14.8 


79.0 


Y(12) 


MWC 1080 


v628 Cas 


23 17 25.59 


+60 50 43.6 


11.6 


7.5 


6.0 


4.8 


15.5+ 1.8 


22.2 


25.1 


Y ( 3) 


Emission-line Stars of Uncertain Classification (possible Herbigs, B 


[e] stars, or other) 












HD 41511 




06 04 59.13 


-16 29 03.9 


5.0 


2.9 


2.1 


1.7 


208.7+20.9 


144.0 


72.2 


N 


HD 45677 




06 28 17.42 


-13 03 11.0 


8.1 


7.2 


6.3 


4.8 


181.9+28.1 


146.0 


143.0 


N 


HD 50138 




06 51 33.41 


-06 57 59.2 


6.6 


5.9 


5.1 


4.1 


94.6+ 9.5 


70.3 


62.5 


N 


MWC 300 


v431 Set 


18 29 25.70 


-06 04 37.2 


10.5 


9.3 


8.2 


6.2 


64.6+ 6.5 


101.0 


113.0 


Y(3) 


MWC 342 


vl972 Cyg 


20 23 03.61 


+39 29 49.9 


10.6 


7.0 


5.8 


4.7 


51.5+13.5 


46.2 


48.6 


N 


T Tauri Stars 
























RY Tau 




04 21 57.40 


+28 26 35.5 


10.2 


7.2 


6.1 


5.4 


14.4+ 1.4 


17.5 


26.1 


N 


T Tau 




04 21 59.43 


+19 32 06.4 


9.6 


7.2 


6.2 


5.3 


10.7+ 1.1 


14.9 


44.2 


Y(13) 


GW Ori 




05 29 08.39 


+11 52 12.6 


9.9 


7.7 


7.1 


6.6 


10.0+ 1.0 


7.9 


20.5 


Y(14) 


FU Orionis Objects 






















FU Ori 




05 45 22.36 


+09 04 12.4 


8.9 


6.5 


5.7 


5.2 


5.9+ 0.6 


5.9 


14.1 


Y(15) 



^Coordinates and JHK magnitudes are from 2MASS (ISkrutskie et al.ll2006h . 
band magnitudes from SIMBAD. 
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^Flux included within aperture of diameter of ?»" and was derived from this work. 

*^IRAS12 and IRAS25 fluxes are from IRAS point source catalog ([joint IRAS Science Working GrouplllQSsh using broadband 
fllters from 8-15/im and 15-30/im. 

^Bi nary references: 1 iKoresko et al.l ([T993|). 2 ICarmona et alJ (l200/t). 3. ICorporon Lagrange! (ll999lV 4. ISmith et aP 
(I2OO5IV 5. iTuthill et all (l200lh. 6. lArava et aD 



Cordero et al ] (l2006l) . 10. [cohen eraD 



Mathieu et al 15. I Wang et al 



7. iThomas et al.l (I2OO7I). 8. Millan-Gabet &: Monnier (2002V 



Monnier et al 



2006 



12. iLeinert et al.. (,19971 13. .Dvck et al.. (4982. ). 



9. 
14. 
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Table 3. 1-D and 2-D Characteristic Sizes at lO.Tytim 



Target 


1-D Gaussian^ 






2-D Gaussian^ 






Comments 




FWHM 


Halo 


xl 


FWHM (mas) 


PA 


Halo 


xl 






mas 


(%) 




Major 


Minor 


deg 


(%) 






Embedded Your 


g Stars (CI 


ass I Sources) 














AFGL 490 


100± 5 


11±3 


7.0 


110± 4 


81± 3 


20± 2 


12±2 


2.4 


Faint extension to South 


Mon R2 IRS 3 


108±16 


39±6 


20.0 


144± 6 


88± 5 


41± 2 


35±2 


6.1 


Nebulosity and binary dX p = 836 mas, 6 = 16.7° 


AFGL 2136 


120± 3 


20±2 


2.1 


125± 5 


115± 3 


44± 14 


19±2 


1.7 


Faint extension to West- Southwest 


AFGL 2591 


118± 3 


9±4 


4.2 


123± 3 


111± 3 


114± 7 


9±4 


3.4 


Nebulosity to West 


S140 IRSl 


134ib 3 


40ibl 


5.5 


136ib 5 


133ib 4 


117ib 25 


40ib2 


5.6 


Extensive nebulosity to South 



Herbig Ae Stars 



AB Aur 


70± 3 


11±1 


0.4 


72± 3 


68± 6 


27± 25 


11±1 


0.4 


Slight extension PA~30° but mc 


MWC 480 


16±20 


3±2 


0.8 


29±24 


0±17 


29± 32 


2±2 


0.8 




HD 142666 


33±24 


4±2 


0.6 


65±65 


0±40 


173± 56 


2±4 


0.5 




HD 142527 


50dz 6 


3±1 


0.4 


55± 8 


45ib 4 


60ib 22 


3dzl 


0.4 




HD 144432 


39dz 5 


Oil 


1.0 


49ibl4 


25ib25 


103± 31 


Oil 


0.9 




MWC 863 


21±21 


2±2 


1.0 


33± 6 


8±14 


12± 21 


2=bl 


0.9 




51 Oph 


0± 4 


Oil 


0.3 


0±27 


0± 1 


93± 53 


Oil 


0.3 




MWC 275 


17±16 


4±3 


1.2 


24±24 


8± 9 


60± 51 


4±3 


1.2 




R CrA 


69± 2 


10±1 


0.9 


77± 2 


64± 2 


78± 4 


10±1 


0.6 




MWC 614 


68± 4 


4±1 


0.7 


80± 3 


52± 8 


17± 9 


4±1 


0.3 




vl295 Aql 


8ibl5 


Oil 


0.5 


34ib34 


0=b 1 


22± 49 


Oil 


0.5 





Herbig Be Stars 



v376 Gas 


37± 3 


1±2 


0.5 


40± 3 


33± 5 


134± 25 


1±2 


0.5 




v892 Tau 


210± 8 


0±3 


21.6 


244± 6 


123± 9 


49± 1 


0±2 


2.1 


Resolved circumbinary disk 


LkHa 101 


70± 3 


4±2 


0.8 


73± 2 


62± 9 


180± 14 


4±2 


0.6 


Slight extension North- South 


v883 Ori 


43± 4 


6±2 


0.9 


46ib 4 


37ibl0 


llOdz 31 


7ib3 


0.9 




MWC 147 


24ibl4 


3ib5 


0.6 


37dzl4 


0±11 


52ib 31 


3ib7 


0.6 




R Mon 


94ib 2 


8±1 


0.7 


97ib 4 


91ib 2 


66± 23 


8=bl 


0.7 


Mostly circularly symmetric 


Z CMa 


61± 4 


1±1 


2.4 


68± 2 


41± 6 


137± 5 


2±1 


0.6 


Faint extension toward companion 


MWC 297 


55± 6 


2±1 


1.2 


59± 9 


50± 6 


40± 21 


2±1 


1.1 




vl685 Cyg 


29±17 


12±6 


0.6 


56±23 


19±19 


112± 26 


11±6 


0.5 




MWC 349 


81± 6 


8±3 


4.0 


94± 7 


26±15 


95± 6 


8±3 


1.1 


Symmetric elongation 


MWC 361 


45±17 


45±2 


0.3 


55±14 


39±30 


17± 26 


45±2 


0.3 


Arcsec-scale halo elongated North- 


v645 Cyg 


46ibll 


5di4 


1.0 


55± 4 


38ibl9 


34± 16 


5±4 


0.9 




LkHa 234 


27±19 


6±1 


0.2 


42±42 


0±44 


91± 34 


6±2 


0.2 




MWC 1080 


0±10 


17±3 


8.1 


0±44 


0± 1 


115± 49 


17±2 


8.4 


Binary: p = 764 mas, = -91.1° 


Emission-line Stars of Uncertain Classification (possible Herbigs, 


B[e] stars, 


or other) 






HD 41511 


Oib 6 


0ib2 


1.6 


19ib26 


9ib 9 


36ib 90 


0ib2 


1.7 




HD 45677 


70ib 3 


12ibl 


1.0 


77dz 3 


66± 3 


70ib 5 


12ibl 


0.7 




HD 50138 


58ib 6 


libl 


1.4 


66ib 4 


46ib 9 


63ib 6 


libl 


1.0 




MWC 300 


49± 3 


1±1 


0.6 


57± 4 


44± 2 


35± 14 


Oil 


0.4 




MWC 342 


38±22 


6±1 


1.0 


39±11 


37±37 


107± 90 


6±2 


1.0 




T Tauri Stars 


RY Tau 


49± 4 


2±1 


0.9 


59±10 


37±25 


12± 45 


2±3 


0.8 




T Tau 


0ibl3 


27ib2 


3.5 


29ib21 


Oib 7 


lib 31 


26±2 


3.6 


Binary: p = 639 mas, = -170.9° 


GW Ori 


37±20 


0ib2 


0.7 


54=b 7 


25±29 


107ib 10 


0ib5 


0.6 





FU Orionis Objects 
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Table 3 — Continued 



Target ' 1-D Gaussian^ 2-D Gaussian^ Comments 





FWHM 
mas 


Halo 
(%) 


xl FWHM (mas) 
Major Minor 


PA 
deg 


Halo xl 
(%) 1 






FU Ori 


21ib22 


10ib5 


1.2 1 55ib27 0ibl3 


65ib 35 


8ib4 1.2 1 Binary: p = 488 mas, = 163.3° 



^Parameters of 1-dimensional Gaussian fits: Full- width at half-maximum (FWHM) in milliarcseconds (mas), percentage of 
light coming from "halo" on scales larger than ~0.5''. 

^Parameters of 2-dimensional Gaussian fits: Full-width at half-maximum (FWHM) along major and minor axes respectively. 
Position angle (PA) of major axis in degrees East of North, percentage of light coming from extended "halo." 



600 r 



pattern6.info 
4x6^IRM OPTIMAL 
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-600 -400 -200 200 400 600 
+X Aperture Plane (cm) 
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15 



10 



> 
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o 

(J) 
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C 
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o 
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O 
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-8 -6 -4 -2 2 4 
Arcseconds on sky (+X) 



Fig. 1. — (left panel) a. This figure shows the 36 segments of the Keck telescope, color-coded by the tilt imposed 
using the Keck Active Control System (ACS) for "pattern 6" described in ^2.11 (right panel) b. This figure shows 
the image plane locations associated with the tilts color-coded in panel (a). The 10^^ X 10^^ array detector can 
simultaneously record fringe patterns of the bottom four image plane locations, while the top image is off the chip. 



— , — I — , — , — , — I — , — , — , — I — , — , 

Segment tilting (Pattern 
Focus-mode Corrected 



■D 
c 
o 
o 



Full Telescope Pupil (focused) 



-2 



-4 



Segment Tilting (Pattern 6) 
Focus-mode Corrected 
+Optical Tweaks 



^ - 




-2 2 

Arcseconds 



Fig. 2. — (left panel) a. This figure shows a single frame (90 ms exposure) of a cahbrator source using the full 
pupil and after careful focusing. Although this image is "diffraction-limited," there are significant variations of this 
PSF with time and seeing conditions and the first Airy ring is usually distorted by phase aberrations, (middle panel) 
b. After applying the tilts and pistons for pattern 6 (see Figured]), we see the light from the star is split into four 
patterns of overlapping fringes. This image was taken after correcting for "focus mode" introduced when changing 
sensor gain setting (see discussion in ^2.2|) . Note the large residual tilt errors on some segments in the top of the 
image frame, (right panel) c. The tilt errors from the middle panel were analyzed at the beginning of the observing 
night and perturbative "optical tweaks" were applied to the segment actuators. The resulting image quality is much 
improved, showing four interference patterns with well-aligned segments. The intensity patterns are all displayed 
with a linear scale. 
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First Measured FWHM Size (mas) 

Fig. 3. — Nearly all the targets of the Keck segment-tilting experiment were observed multiply times under in- 
dependent observing conditions - this was a core design principle of our observing methodology. Here we plot the 
measured size (Gaussian FWHM) at one epoch versus another in order to estimate systematic errors as a function of 
the target size. Here we see extremely good self-consistency within reported errors and the typical 1-sigma confidence 
limits for the mean diameter are shown in the central grey-shaded region. We can reliably (with 2-sigma confidence) 
estimate partially-resolved objects down to a full- width at half- maximum of ~35 milliarcseconds, about 10 x smaller 
than the formal diffraction limit. 



CM 



■D 


CO 

c/D 



"co 
> 



1.5 
1.0h 

0.5- 

0.0 
1.5 

1.0 

0.5 

0.0 
1.5 

1.0 

0.5 ^ 
0.0 



Embedded Young Stars: 
AFGL 490 





®s ■- _ J 



halo:10% 
fwhm:100 mas 



Herbig Ae Stars: 
AB Aur 



Herbig Ae Stars: 
MWC 863 



halo: 2% 
fwhm: 20 mas 



Embedded Young Stars: 
Mon R2 IRS 3 




Herbig Ae Stars: 
MWC 480 



Herbig Ae Stars: 
51 Qph 



halo: 0% 
fwhm: mas 



Embedded Young Stars: 
AFGL 21 36 



halo: 19% 
fwhm:1 19 mas 



Herbic) Ae S;rs;: 
HE)142()f)6 




Herbig Ae Stars: 
MWC 275 




halo: 4% 
fwhm: 16 mas 



Embedded Young Stars: 
AFGL 2591 




halo: 8% 
fwhm:1 18 mas 



Herbig Ae Stars: 
HQ- 142527 



Herbig Ae Stars: 
ROrA 



halo: 10% 
fwhm: 69 mas 



Embedded Young Stars: 
S140 IRS1 



halo:39% 
fwhm:134 mas 
-h 



Herbic Ae Stars: 
HQ)1 44432 



Herbig Ae Stars: 
MWC 614 



halo: 4% 
fwhm: 67 mas 



0246802468024680246802468 

Baseline (m) 



1.5 
1.0 

0.5 

0.0 
1.5 

1.0 

0.5 

0.0 
1.5 

1.0 

0.5 
0.0 



Herbig Ae S^^r^: 
vt2J 



Herbig Be Stars: 
v376 Cas 



Herbig Be Stars: 
v892 Tau 




halo: 3% 
fwhm: 24 mas 



halo: 7% 
fwhm: 93 mas 



halo: 1% 
fwhm: 60 mas 



Herbig Be Stars: 
LkHa101 




halo: 3% 
fwhm: 69 mas 



Herbig Be Stars: 
MWC 297 



C) 




halo: 1% 
fwhm: 54 mas 




Herbig Be Stars: 
v883 Ori 





Herbig Be Stars: 
MWC 349 




Herbig Be Stars: 
MWC 361 



Herbig Be Stars: 
v645 Cyg 



Herbig Be Stars: 
MWC 1080 



halo: 8% 
fwhm: 80 mas 






halo: 5% 
fwhm: 45 mas 



0246802468024680246802468 

Baseline (m) 



"D 
CD 



> 




2 4 6 8 



024680246802468 

Baseline (m) 



Hi 

o 
El 

CD 
bO 



0.0 0.2 0.4 0.6 0.8 1.0 



Embedded Voung Stars: 



AFGL 490:: Mon R2 IRS 3 



Herbig Ae Stars: 



Embedded Voung Stars: 



Herbig Ae Stars: 
51 Oph 



Embedded Voung Stars: 
AFGL2136i 



Visibility-Squared 




Herbig Ae Stars: 



mbedded Voung Stars 
AFGL 2591 1 



Werb'ig Ae Stars: : 
RCrA 



Embedded Voung Stars: 
S140 IRS1 



Werb'ig Ae Stars: 
614 



'big Ae ; 
MWC 



863 



■7 




10 5 -5 -1010 5 -5 -1010 5 -5 -1010 5 -5 -1010 5 -5 -10 

u coordinate (meters) 



0.0 



0.2 0.4 



0.6 0.8 1.0 



I ' ' ' ' I ' 



Herbig Ae Stars: 
V1295 Aql 



Herbig Be Stars: 
MWC 147 



I ' ' ' ' I 



Herbig Be Stars: 
v376 Cas 



Herbig Be Stars: 



Visibility-Squared 



Herbig Be Stars: 
v892 Tau 



Herbig Be Stars': 



Herbig Be Stars 
LkHa101i 




Herbig Be stars': 



— p-i— I— I— I— p-i 



Herbig Be Stars: 
v883 Ori 



Herbig Be Stars: 
V1685 Cyg 



R Mon 



ZCMa 




IHerbig Be Stars: 
349 



Herbig Be Stars: 





Herbig Be Stars: 
v645 Cyg 



Herbig Be Stars: : 
LkHa 234: 



Herbig Be Stars: 
MWC 1080 








10 5 -5 -1010 5 -5 -1010 5 -5 -1010 5 -5 -1010 5 -5 -10 

u coordinate (meters) 



Uncertain Classificatioh: 
HD41511 




Uncertain Classificatioh: 
HD 45677 




Visibility-Squared 



Uncertain Classification: 
HD 50138 




Uncertain Classification: 
MWC 300 




'F^ubrionis'6bi iects: 



'Uncertain Classiflcatioh: 
MWC 342 




a; 

Hi 

o 

CD 
bO 



"Tauri Stars: 
RY Tau 



"Tauri Stars: 
TTau 



" Tauri Stars: 
GW Ori 







10 5 -5 -10 



10 5 -5 -10 



10 5 -5 -10 
u 



10 5 -5 -1010 5 -5 -10 
coordinate (meters) 



-35- 



0.010 0.032 



0.10 



0.32 



1.0 



Relative Intensity 




■ Mon R2 IRS 3 



0.5 0.0 -0.5 



Mon R2 IRS 3 
(zoomed) 



AB Aur 
(zoomed) 




0.3 0.2 0.1 0.0-0.1-0.2-0.3 



0.3 0.2 0.1 0.0-0.1-0.2-0.3 



0.3 
0.2 
0.1 
0.0 
-0.1 
-0.2 
-0.3 



ZCMa 
(zoomed) 



0.3 
0.2 
0.1 
0.0 
-0.1 
-0.2 
-0.3 



MWC 349 
(zoomed) 



1.0 
0.5 
0.0 
-0.5 
-1.0 



MWC 1080 



0.3 
0.2 


AFGL2136 
(zoomed) 


0.3 
0.2 


AFGL2591 
.(zoomed) 


0.1 




0.1 




0.0 




0.0 




-0.1 




-0.1 


■ 


-0.2 




-0.2 




-u.o 




-0.3 




0.3 0.2 0.1 0.0-0.1-0.2-0.3 


0.3 0.2 0.1 0.0-0.1-0.2-0.3 


0.3 
0.2 


LkHalOl 
.(zoomed) 


0.3 
0.2 


R~Mon 


0.1 
0.0 
-0.1 


# : 


0.1 
0.0 
-0.1 


w 


-0.2 




-0.2 




-0.3 




-0.3 




5 


0.3 0.2 0.1 0.0-0.1-0.2-0.3 


0.3 0.2 0.1 0.0-0.1-0.2-0.3 


1.0 


TTau 


1.0 


FU Ori 


0.5 




0.5 




0.0 


• 


0.0 




-0.5 




-0.5 




-1.0 




-1.0 





0.3 0.2 0.1 0.0-0.1-0.2-0.3 



0.3 0.2 0.1 0.0-0.1-0.2-0.3 1.0 0.5 0.0 -0.5 -1.0 1.0 0.5 0.0 -0.5 -1.0 1.0 0.5 0.0 -0.5 -1.0 

Arcseconds 



Fig. 6. — This figure shows the aperture synthesis images of the most resolved targets in our sample. The intensity 
scale is logarithmic and is specified in the color bar. North is up and East is left in all figures. 
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Fig. 7. — This figure shows imaging of MWC 361 A at 1 1.6/im and 17.65/im using the LWS camera on the Keck 
telescope without segment-tilting (first presented in Perrinl t2QQ6). We see an arcseco nd- scale North- South extension 
which follows the orientation of the underlying sub- arcsecond binary in this system ( Monnier et al. 2006h . This is a 
remarkable structure suggesting the circumbinary disk is in an advanced stage of photoevaporation. The intensity 
scale is logarithmic and is specified in the color bar. North is up and East is left in all figures. 
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Stellar Luminosity (solar) 

Fig. 8. — The mid- infrared size - luminosity diagram. Different plot symbols correspond to different YSO types: 
Embedded YSOs (asterisk), Herbig Ae (X), Herbig Be (diamond), T Tauri (triangle), FU Orionis objects (square), 
and non- classified emission line objects (plus). System s in close binaries are labeled (b). The lines represent different 
temperatures for grey dust according to definition in iMonnier Millan-Gabet (2002 ). Objects with FWHM sizes 
smaller than 35 mas are labeled as upper limits. 



- 38 - 



Angular Size vs IRAS Color 



9K Embedded 
X Herbig Ae 
Herbig Be 
+ Uncertain 
A TTauri 
□ FU Ori 




Luminosity-normalized Ring Radius vs IRAS Color 



Embedded 
X Herbig Ae 

Herbig Be 
+ Uncertain 
A TTauri 

FU Ori 



I ,T, 



I I 



Fig. 9. — (left panel) a. The mid- infrared angular size (Gaussian FWHM) vs. IRAS 12-25 color (-2.5 log F^(12/im) 
/ Fiy(25/im) ). Systems in close binaries are labeled (b). Objects with FWHM sizes smaller than 35 mas are shown 
as 35 mas upper limits, (right panel) b. Here we plot the luminosity- normalized ring radius (AU/root luminosity) 
as a function of the IRAS 12-25 micron. The luminosity-normalized ring radius is independent of distance and is 
diagnostic of disk structure. Note the significant correlation for the Herbig Ae sample. 



